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ABSTRACT 

The spin of Cygnus X- 1 is measured by fitting reflection models to Suzaku data covering the 
energy band 0.9-400 keV. The inner radius of the accretion disc is found to lie within 2 grav- 
itational radii (r g = GM/c 2 ) and a value for the dimensionless black hole spin is obtained of 
0.97^o Q 2 4 . This agrees with recent measurements using the continuum fitting method by Gou 
et al. and of the broad iron line by Duro et al. The disc inclination is measured at 23.7^‘J deg, 
which is consistent with the recent optical measurement of the binary system inclination by 
Orosz et al of 27 ±0.8 deg. We pay special attention to the emissivity profile caused by irradi- 
ation of the inner disc by the hard power-law source. The X-ray observations and simulations 
show that the index q of that profile deviates from the commonly used, Newtonian, value of 3 
within 3 r g , steepening considerably within 2r g , as expected in the strong gravity regime. 
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1 INTRODUCTION 

Astrophysical black holes are characterised by just mass and spin. 
The measurement of spin requires observations which probe the 
immediate environment of the black hole event horizon. Two meth- 
ods which have emerged using cuirent X-ray observations of lu- 
minous accreting black holes rely on fitting either the shape of 
the reflection component of the spectrum including a broad iron 
line, or the shape and the flux of the quasi-blackbody continuum. 
Both assume that the accretion disc extends in to the Innermost 
Stable Circular Orbit (1SCO), the radius of which is determined 
by the black hole spin. The Reflection Fitting method essentially 
measures the largest gravitational redshift of the disc, whicb comes 
from the ISCO, and yields that radius in units of gravitational radii 
(r g = CM /(?■). The Continuum Fitting method measures the area 
of the inner disc and obtains the radius of the ISCO in km. Here we 
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apply the reflection method to Suzaku X-ray data of the first stellar 
mass black hole, Cyg X-l. 

Early attempts to measure the spin parameter of Cyg X-l us- 
ing the reflection method led to inconsistent results (Miller at al 
2005; Miller et al 2009). Adopting a Newtonian emissivity profile 
(surface flux on the disc varies with radius r as r~ q with q = 3) 
and using a special timing mode of XMM, a recent measurement 
by Duro et al. (201 1) gives a dimensionless spin parameter 1 of 
a — 0.88io‘n- results using the Continuum Fitting method 

have been obtained by Gou et al (2011). They find a near ex- 
treme black hole with a > 0.95. The work relies on accurate op- 
tical determinations of the mass (14.8 ±0. 1 M©) and inclination 
(27 ±0.8 deg) obtained by Orosz et al (201 1). 

Here we study the reflection/iron line approach in detail, tak- 
ing account of the relativistic effects on the emissivity profile ex- 


1 The dimensionless spin parameter a — cJ/GM 2 , where J and M are the 

angular momentum and mass of the black hole, respectively. 
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pected if the black hole spins rapidly. Such effects demand that 
the emissivity be steep if the disc extends to the innermost region 
within ~ 2 r g (Wilkins & Fabian 2011). We find a = 0.97t[};Q2 4 * 
which is a higher than, but marginally consistent with, that reported 
by Dura et al (201 1). Both reflection results are consistent with the 
Continuum Fitting result of Gou et al (2011). 

We find the disc to be fairly highly ionized with an ionization 
parameter around 1800 erg cm s" 1 , which means that the reflec- 
tion spectrum is dominated by a large edge due to ionized iron. 
This blurred edge carries information about spin in the spectral fit. 
Early work emphasised how this edge can be important in fitting 
the spectrum of Cyg X-l (Ross, Fabian, Brandt 1996). 


2 SPECTRAL FITS 

The dataset we use is one of the 20 Suzaku observations of Cyg 
X-l analyzed by Miller et al (2012; the observation of 2009 Apr 
8). The focus of that work is not on measuring the black hole spin, 
but rather studying the detailed disc-jet coupling in Cyg X-l. The 
spectra in that work have been fitted with a relativistically-blurred 
reflection model, but in order to make the work tractable, the fits 
kept the source inclination and iron abundance frozen at 27 deg 
and 1, respectively. The emissivity profile of the blurring function 
was also just a single power-law, fixed for the Newtonian value of 
3. We relax these assumptions here and apply them to an average 
dataset. We later tested the model on several other datasets (2009 
May 6, 25 and Dec 1) and obtained consistent results. 

The best fit to the 1-500 keV spectrum 
is shown in Fig. 1 and values of key parame- 
ters are given in Table 1. The model used, CON- 
STANT -TBABS *MTABLE(WINDABS. FITS) (GAUSSIAN+GAUSSIAN 

+diskpbb+powerlaw+kdblur2f : atable(extendx.mod)) 
highecut, is shown in Fig. 2. KDBLUR2F is a fast adaptation 
by JSS of the broken power-law emissivity convolver kdblur2. 
Extendx is a version of the self-consistent slab reflection model 
REFLIONX (Ross & Fabian 2005) extended to allow the photon 
index to be less than 1.4. (We do not use the REFHIDEN models of 
Ross & Fabian 2007, which incorporate a blackbody component 
because they are only available for Solar abundance iron and at 
the blackbody temperature found here for the low state of about 
0. 16 keV the differences are minimal.) Absorption features due to 
the stellar wind of the companion (e.g., Hanke et al. 2009, Nowak 
et al., 201 1) are modeled using windabs.fits (Miller et al., 2012) 
and we refer to that paper for further details. The absorption has 
little effect on our results. We note that the highecut model 
used to fit the data above 100 keV is phenomenological and is not 
consistent with the 300 keV exponential cutoff assumed by the 
EXTENDX model, but has limited effect on the results presented 
here. The quality of the fit is very good, and the emissivity index 
q is found to be high > 6.8. The inclination of 23.7+S?J deg is 
consistent with the optical measurement of Orosz et al (2011). 

It is clear that the spectrum is reflection-dominated below 
10 keV (Fig. 2), although the whole source, from an energy point 
of view integrating over the whole spectrum, is not reflection dom- 
inated. The ratio of the 0.1-1000 keV flux from the reflection com- 
ponent is 0.70 times that from the power-law component, so the 
reflection fraction 3$ = 0.7. The blackbody disc component is 0.4 
times the power-law and if reflection is responsible for some of 
that emission then 3Z ~ l. A narrow emission line is introduced 
at 6.4 keV to account for (weak) distant reflection (e.g. from the 
stellar companion) and a narrow absorption feature is included at 



Energy (keV) 

Figure 1, Spectrum of Cyg X-l in the Suzaku XISOt-XIS3 (black), PIN 
(red) and GSO (green) detectors, from Miller et al (2012), The ratio to the 
best fitting model (Table 1) is shown below. The data have been rebinned 
for display purposes only. 



Figure 2. Best-fitting spectral model plotted as EFe. The relativistically- 
blurred reflection component is shown as the red curve. The 100 keV kink 
introduced by HIGHECUT has no effect on our spin results. 


1 .85 keV to compensate for problems in the model response near 
the detector silicon edge. 

We show the model spectrum over the 2-65 keV band in the 
top panel of Fig. 3. The inner radius has been stepped from 1.3 to 
1.9, and then to 400r g and the innermost emissivity index q — 9. 
The effect of fixing q at 3 is shown in the lower panel: the model is 
now less sensitive to the innermost tadii. 

We next fit the spectrum with a model in which the blurring is 
carried out on many small contiguous annuli, simultaneously. This 
is similar to the approach taken by Wilkins & Fabian (201 1) for the 
AGN 1H 0707-495. Once again, the emissivity has a steep index at 
small radii before breaking to a flatter part, then dropping further 
close to the Newtonian expected value of 3. Simulations produced 
by ray tracing in the strong gravity regime revealed in Wilkins & 
Fabian (2012) show that this triple power-law shape is expected 
from an irradiated disc around a rapidly spinning black hole. The 
outer break suggests that the height of the irradiating source above 
the disc is around 5 — 7r g . The results show that emission from 
within 2r g is definitely required. 

The model which was used to fit the data from the XIS, 
PIN and GSO simultaneously uses kdblur2f, which is a broken 
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Figure 3. (Top) Model reflection spectra, relativistically blurred for inner 
radii of 1 ,3, 1 .5, 1 ,7, 1 .9 and 400 r g with an emissivity index q of 9. (Bottom) 
Similar to above but with q = 3 and expanded scales. Note the reduced 
sensitivity to small values of r m . 



Figure 4. Emissivity profile obtained by fitting the blurred reflection model 
to small increments in radius. The profile has the characteristic steep - flat- 
moderate profile expected from irradiation by a source at a height h-5- 
7r g above die black hole. The innermost annulus with a significant detection 
is at r ~ 1.4r g . 


Thble 1 . Values of key model parameters used in the spectral fitting. 


Component 

Parameter 

Value 

q fixed 

powerlaw 

Photon index, T 


1 *^- 0.009 

kdblur 

Inclination, deg 
Rio. Tg 

Index, qi 
Index, <72 

break* r g 

23.7 IK 

A ' W - 0.16 

>6.8 

2.75 ±0.15 
4.0±1.1 

39.8^2 

<1.66 

3 

3 

reflionx 

Iron abundance / solar 
Ionisation parameter, £ 
ergcms“ ! 

1 40 + 0.64 

1 — 178 

0.89*0.12 

29601“ 


X 2 / dof 

2217/2419 

2411/2421 



Figure 5. Effect of varying the inner radius as a function of iron abun- 
dance Fe. 


power-law capturing the essence of Fig. 4, We investigated the de- 
pendence of the inner radius on various key parameters* the iron 
abundance and the inclination (Figs. 5-7). 

The data require an iron abundance which is slightly superso- 
lar, an inclination less than ~ 30 deg and an intermediate ionization 
parameter around 1760. All results strongly point to an inner radius 
less than 2 r gi even if q is fixed at 3 (Table 1, Fig. 7). 

As a test, we have fitted the PIN data alone, since this dataset 
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is sensitive to the reflection hump. The iron abundance is fixed at 
unity (the depth of the edge in the XIS data indicates that the abun- 
dance cannot be either very low or very high). We also fix the in- 
clination at 27 deg, obtaining again r in < 2r g . This confirms that 
fitting the reflection hump alone can in principle measure spin. 

The reflection fits presented here all strongly point to the spin 
of the black hole being high. Using the broken power-law emis- 
sivity model model relconvf (a faster modification by JSS of 
RELCONV of Dauser 2010, 2012) yields a — 0.97^q ^ 4 at the 90 
per cent confidence level. The results for the other parameters are 
similar to those shown in Table 1. 

We have also used the Novikov & Thome (1973; see also 
Page & Thome 1974) emissivity profile which is appropriate if 
the corona is immediately above the disc. The emission then varies 
with radius in the manner expected from a radiatively-efficient thin 
disc. The resulting fit has x 2 = 2422/2422, significantly worse than 
the earlier best fit. High spin (a > 0.98) is again required and the 
inclination is higher (/ ~ 39 deg). 


2.1 The Newtonian power-law emissivity profile, q — 3 

Duro et al (201 1) measure the broad iron line of Cygnus X-l as seen 
by XMM-Newton in the EPIC-pn modified timing mode. They re- 
port that the black hole is spinning close to maximal with value 
a — 0.881q jp on adopting the Newtonian, single power-law emis- 
sivity profile with index fixed at q = 3. A fit with q free gave 
q = 10ig and a — -0. 1 ± 0.4. They rightly dismiss the statistically 
more probable, high q - low spin, solution as making no physical 
sense, since high q is only applicable at high spin. 

Fitting our spectra with q fixed at 3 leads to a much worse 
(although acceptable) fit with A# 2 — 4- 194 (Table 1). The inner ra- 
dius is inferred to be at less than 1.66r g , the inclination is above 
35 deg (which is inconsistent with the optical result) and the ion- 
ization parameter is higher at nearly 3000. If the inclination is fixed 
ar 25 deg, then r m < 1.9r g . The conclusion that the black hole is 
rapidly spinning remains robust. 

The spectral residuals obtained from using a single fixed value 
of q and a broken power-law with variable indices are shown in 
Fig. 7. The large difference in x 2 is mostly due to the red wing of 
the line feature at 2.5-3 A keV and (to a much lesser extent) the 
edge at 6.5 — 7.5 keV. The residuals for fixed q are at about 2 per 
cent and increase to 3 per cent or more if r, n is increased to 3r g and 
beyond. This demonstrates where our information is coming from 
and shows that from a statistical point of view we can probe the 
strong gravity regime. 

A single power-law emissivity has only limited validity. It is 
clear from the fitted emissivity profile in Fig. 4 that a slope of 3 
is a fair approximation from r = 2r g outward. It underestimates 
the profile within r — 2r g and is therefore a poorer probe of the 
innermost region (see Fig. 3). If q = 3 is used, then it will likely 
yield an upper limit to the inner radius and thus a lower limit on the 
spin. Doubly broken power-law profiles are required when fitting 
for the inner radius and/or spin when the spin is high (r ra < 2r g and 
a > 0.94). 

Cyg X-l is a good example where a steep inner profile is 
needed (Fig. 8). In contrast, the source GX339-4, which has an in- 
ner radius of ~ 2r g and thus moderate spin, requires only a mildly 
steep index q of around 4 (Fig. 7). (We use the XMM data of the 
low-state observation of GX339-4 from Reis et al 2008 for this fig- 
ure.) 

The higher state of Cyg X-l during the Duro et al (2011) ob- 



Figure 7. Ratio of (heavily binned) spectral daiapoints to the best fit- 
ting model with a single power-law emissivity profile (index q fixed at 3: 
black points) and to the broken power- law model with variable indices (red 
points). 

servations, in which the black body component contributes to about 
3 keV, may complicate variable q measurements. 


3 MODIFICATIONS TO ALLOW FOR THE TIME 
VARIABILITY OF CYG X-l 

In this Section we consider what modifications can be made to the 
spectral model developed here to allow for the complex time vari- 
ability observed in Cyg X-l. The aim is not to build a complete 
model which accounts for everything that has been observed but to 
indicate the directions in which modifications are needed. 

Study of the rapid variability of Cyg X-l requires high count 
rates and thus instruments with large collecting areas, such as EX- 
OS AT (Belloni & Hasinger 1990), Ginga (Miyamoto et al 1988, 
1989; Gicrlinski et al 1997) and RXTE (Nowak et al 1 999; Revnivt- 
sev, Gilfanov & Churazov 1999; Gilfanov, Churazov & Revnivtsev 
2000; Pottschmidt et al 2003). The power spectrum is flat up to 
about 0.02 Hz, above which it drops with frequency / as / * up to 
a few Hz where it steepens to f~ 2 . Time lags are seen over the fre- 
quency range 0. 1-30 Hz in which soft bands lead hard bands, with 
the longest time lags being 0.05 s. As a light crossing time, this 
corresponds to a distance of ~ 1000r g , which is much larger than 
the radius of any of the regions inferred by our spectral fitting. The 
lag timescale must therefore be some form of propagation time. In 
the low state, such as observed here, the iron “line’* is found to be 
variable up to a frequency of a few Hz (Revnivtsev, Gilfanov & 
Churazov 1999). This is too low a frequency to be compatible with 
any light crossing time effects expected in our spectral model. 

Several models have been made for the spectral variability and 
time lags. Some early models are reviewed by Poutanen (2001). 
A model by Lyubarskii (1997) in which fluctuations are generated 
by and propagate inward through the disc has been extended by 
Arevalo A Uttley (2007) to explain many of the features observed. 
Most recently the hard state lags have been modelled as propagat- 
ing oscillations in a hot inner flow in the central region of an accre- 
tion disc truncated at about 20r g (Kawabata & Mineshige 2010), 
prompted by results of an earlier Suzaku observation of Cyg X-l 
reported by Makishima et al (2008). In that work the corona was 
modelled as having two optical depths, which allow for softer and 
harder continue. Reflection was dealt with by a phenomenological 
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Figure 8. Variation of inner emissivity index q\ on inner radius r m for Cyg 
X-l (tcp) and GX339-4 (bottom). A broken power-law emissivity profile is 
used with the outer index, q 2 , fixed at 3, for both plots. 



Figure 9. Emissivity - spin contours for (left, black) a simple power-law 
profile and (right, red) a broken power- law profile with the outer index fixed 
at 3. The minimum x 2 for the fit using the broken power-law is about 150 
be ow that for the fit using a single power-law emissivity profile. The results 
in Tkble 1 were obtained with the outer index free. 


Ihble 2. Values of model parameters obtained from the best fit with two 
power-law components of spectral indices T s and iy The incident spectrum 
for the reflection has index TV 


Component 

Parameter 

Value 

powerlaw 

Photon index, r 8 

1.45f 

powerlaw 

Photon index, Th 

1.25f 



1.63 

reflionx 

Iron abundance / solar 

1.28 


Ionisation parameter, £ 
erg cm s^ 1 

1775 


* 2 /dof 

2219/2419 



Figure 10. Model spectrum with two power-law components. The index of 
the incident spectrum for the reflection component equals that of the steeper 
power-law. 


model involving cold unblurred reflection and abroad gaussian iron 
line. 

The results obstained from the fits presented here reveal an in- 
ner radius of less than 2r g strongly irradiated by a corona situated 
at ~ 5 — 7 r g . This is incompatible with an extensive inner corona. 
The single power-law continuum used here does not allow for time 
lags unless it varies in slope on the required timescales. To allow 
for slope variations we have fitted the spectrum with two power-law 
continua, the steeper of which gives rise to the reflection. The spec- 
tral fitting is degenerate provided that the sum of the two power- 
laws approximates the single power-law shown in Fig. 2. For illus- 
tration purposes we have selected a pair of spectral indices which 
give a good fit (Fig. 10; Table 2). Reflection continues to dominate 
the spectrum below 10 keV. 

A simple physical interpretation arises if the source of the 
harder power-law is quasi-static and situated close to the centre of 
the disc, therefore producing the reflection. The steeper component 
could originate from further out within the moving jet; any reflec- 
tion produced by this continuum is reduced due to beaming (Be- 
loborodov 1999). Time lags then occur if variations propagate up 
the jet from the inner, quasistatic and softer component to the outer, 
jetted harder component. They need not propagate at jet speed, but 
could be due to changes in the magnetic structure or collimation of 
of the jet. The power in these components and the jet is presumed 
to derive from the accretion disc and be transferred by magnetic 
fields coupled to a range of radii, so varying in a manner which is 
consistent with the Lyubarskii model. 
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The lack of rapid reflection variations reported by Revnivtsev 
et al 1999) could in part be due to light bending effects implicit 
in the proximity of the inner, softer, power-law component to the 
black hole (Miniutti & Fabian 2004), if some of the rapid variability 
is due to (small) changes in height of that inner component. 

Further development of such models is beyond the scope of 
this paper, which is focussed on the low state X-ray spectrum of 
CygX-1. 


4 DISCUSSION 

We have used the X-ray Reflection Fitting Method to obtain a ro- 
bust measurement of the dimensionless spin parameter of the black 
hole in Cyg X-l of a — 0.97l^ ^ 4 . Our result agrees with the con- 
tinuum fitting value of Gou et i (201 1) and yields an even higher 
spin than the reflection measurement of Duro et al (2011). Cygnus 
X-l definitely appears to have high spin. 

We argue that non-Newtonian values must be used for the in- 
ner emissivity profile of the reflected emission of the innermost disc 
when the spin is high (r m < 2r g and a > 0.94). 

The Reflection Fitting method relies on detection of the in- 
nermost radius of the disc which is dense enough for reflection to 
occur. We assume that this is the ISCO, for within that radius the 
plunge orbits mean that the gas density drops rapidly so the reflec- 
tion signal vanishes. Although magnetic fields may have some ef- 
fect here, we consider it reasonable to assume that the dense parts 
of the disc that give the reflection spectrum are not seriously af- 
fected (see discussion and simulations in Reynolds & Fabian 2008 
and Shafee et al 2008). There is nevertheless a systematic uncer- 
tainty here which simulations need to resolve. 

The reflection fraction & ~ 0.7 is low compared with the sim- 
ple expectation of a source close to the black hole at h ~ 5r g . Light 
bending should give ^^2 (e.g. Fabian et al 201 1). The observed 
radio emission (Miller et al 2012) means that a jet is operating, 
so it is plausible that the emission model is more complex, as dis- 
cussed in Section 3. We suggest there that irradiation of the disc 
by the quasi-static base of the jet at a 5 - 7r g gravitational radii is 
responsible for the reflected emission, whereas the observed power- 
law continuum is dominated by higher, faster parts of the jet which 
mildly beam emission in our direction. 

Systematic uncertainties in the models and (relative) spectral 
calibration can also affect the details of our fits and thus results, par- 
ticularly of the innermost region. The Continuum Fitting method is 
also subject to systematic uncertainties in both models and spectral 
calibration (relative and absolute). 

Wide agreement on the high spin of Cygnus X-l boosts con- 
fidence in spin measurements from the Reflection Fitting method 
of other objects as well as AGN, where it is the only method yet 
available. 

The spectral fits strongly indicate that the innermost radius of 
the disc lies within 2r g . Relativistic effects are strong there and a 
simple Newtonian assumption for the emissivity profile is inappro- 
priate. The effects are inescapable and involve both special relativ- 
ity, since the velocity of the reflecting matter on the disc is high, and 
general relativity, through gravitational redshift and light bending 
(Wilkins & Fabian 2010, 2012). 

Theoretical emissivity profiles due to point sources at differ- 
ent heights on the rotation axis above the disc plane are shown in 
Fig. 12 of Fabian et al (201 1). When the source height exceeds 3r g , 
the emissivity profile breaks at radius r ~ h to the Newtonian value 


of 3 outside and flatter within. The inner profile rapidly steepens to 
a high value within 2r g , if the disc extends in that far. 

We recommend that q = 3 be used initially. Should this lead 
to high spin being suspected then a power-law broken at least once, 
with a high inner index (q\ > 3), should be used. This is illustrated 
in Fig. 9 where we have fitted our Cyg X-l data by a single power- 
law emissivity of free index. It shows the best fit being of high index 
and relatively low spin ( q ^ 10, a ~ 0.6), which is rejected by the 
physical sense argument of Duro et al (2011). If instead we use a 
broken power-law emissivity fixing the outer index qi = 3, then we 
find that the inner emissivity q\ remains high and the spin rises to be 
very close to maximal ( a > 0.95). This now makes physical sense. 
Allowing the outer index to be free gives the results of Table 1. 
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